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This constitutes the semi-annual status report for the period 
July 1, 1982 - Dec. 21, 1982, on the research being performed by the 
School of Aeronautics and Astronautics, Purdue University, for the NASA 
Dryden Flight Research Facility, Ames Research Center, under Grant NAG4-1. 
The objective of this research effort has been the development of a uni- 
fied control synthesis methodology for complex and/or non-conventional 
flight vehicles, and to understand, enhance, and develop prediction techn- 
niques for the handling characteristics of such vehicles. 

2. Publications, Personnel and Discussion 

Two papers were presented at the 1982 AIAA Atmospheric Flight Mecha- 
nics, and Guidance and Control Conferences in San Diego in August, 1982, 
on results obtained under this grant. These papers were 

1. "Application of An Optimal Cooperative Control Technique for 
Augmentation Synthesis of a Control Configured Aircraft," by 
Mario Innocenti and David. K. Schmidt, presented at the Guidance 
and Control Conference AIAA Paper No. 82-1520. 

2. "A Modern Approach to Pilot/Vehicle Analysis and the Neal-Smith 
Criteria," by Barton J. Bacon and David K. Schmidt, presented 
at the Atmospheric Flight Mechanics Conf. , AIAA Paper No. 82- 
1357. 

In addition, the material documented in NASA CR 163112, entitled 
Pilot-Orientation Multivariable Control Synthesis by Output Feedback was 
submitted over a year ago for possible publication in the AIAA Journal 
of Guidance, Control and Dynamics. The reviewers comments have been 
received during this reporting period, with their recommendation that the 
paper be accepted for publication after the material in paper 1) above is 
added to the manuscript. This is currently being accomplished and we are 
eager to have these results appear in a journal article. 
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Also, paper 2) above was submitted as well for journal publication. 
The reviewers have recommended that this paper appear after only minor 
revision in the Journal of Guidance, Control and Dynamics. This has 
been accomplished, and we look forward to its appearance. 

Finally, a paper has just been completed, and will be presented at 
the 1983 IEEE National Aerospace and Electronics Conference (NAECON), 
Dayton, OH, May 1983. This paper is entitled 

"Integrated Pilot-Optimal Augmentation Synthesis for Complex 

Flight Systems: Experimental Validation," by Mario Innocenti 

and David K. Schmidt. 

This last paper documents the experimental validation of the analytical 
predictions of the handling characteristics of several sets of augmented 
vehicle dynamics. Copies of tho manuscripts for all the papers cited 
above have been previously forwarded to the Technical Monitor. 

With the graduation of Mario Innocenti in December, 1983, extensive 
development of the cooperative control synthesis theory is complete, and 
future activities will most likely focus on application of the methodo- 
logy to a variety of areas. 

The research on pilot model identification and on extending the 
Neal -Smith approach is progressing well. Mr. Pin- Jar Yuan is developing 
several approaches to pilot model identification, extending the work in 
this area that has been reported on in several previous papers. The 
status on this work is completely documented in the Appendix to this 
report. 

The investigation of the potential of extending the optimal -control 
approach to Neal/Smith analysis for the approach and landing task is the 
final area of research currently being pursued. This work is being per- 
formed by Mr. Bart Bacon, the co-author of our paper on the method, and 
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Mr. Dan Garrett, a new student now supported on this grant. Mr. Garrett 
is an M.S. student and is planning to graduate in December, 1983. This 
topic area will constitute his M.S. thesis. 
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b Background 
!• b Introduction 

* 

In the late 60's# the optimal control model (OCM) of 
the pilot . was developed* uihich is based on the hypothesis 
that a human operator (pilot) chooses his control input to 
minimize some cost function subject to his known physical 
limitation. The OCM of the human operator has yielded 
results for the manual control of a variety of plants that 
agree with experimental findings provided that the correct 
cost function is assumed (Ref. 1). 

Later* Hess (Ref. 2) shouted that there exists a strong 
correlation between the subjective pilot evaluation (e. g. 
Cooper -Harper rating) of the vehicle and the magnitude of 
the OCM quadratic cost function. Recent investigations 
(Ref. 3) have provided more substance to the idea that such 
a correlation exists over a wider variety of piloted vehicle 
tasks. The above correlation between the pilot rating and 
the OCM objective function has been used by Schmidt (Ref. 
4*5) in the attempt to develop a unified theory of vehicle 
handling qualities and optimal flight control synthesis. 

JL 2, Model structure 

The analysis relies on the well-known (Ref. 1) optimal 
control theoretic technique for modeling the human pilot 
manual control function. The hypothesis upon which it is 
based is that the well trained* well motivated pilot chooses 
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OF 


hi* control input* (e. g. stick force) to meet the pTTot's 
mission objective/ which can be described a* minimizingg a 
meet the pilot's mission objective# which can be described 
as minimizing a suitable cost function in the task# subject 
to his human, limitations. The cost function i* further 
assumed to be expressible in terms of a quadratic form as 
following; 


j t | 

J « E I lim J- X < y T ay + u T Ru > dt I 

I t~>oo 0 I 


<1 . 1 > 


where 


y = vector of pilot's observed variables 
u « vector of pilot's control inputs 
Q# R - pilot selected weightings 

and the pilot model is sketched briefly in Figure 1-1. The 
human pilot chooses his "best” control decision (for example 
stick force) based on the information displayed to him 
(pilot observations) and the performance objective refer- 
ence. So the suitable selection of cost function is very 
important in representing the human pilot. And in this pilot 
model# we include the human limitation such as information- 
acquision* time delay« observation and control output 
noises* and neuromuscular dynamics etc. . 

The pilot perceives measurable variables y delayed a 

fixed time r , and contaminated by white observation noise 

v . So his measurement vector u is 
y 


y p <t> = y<t-r> + v y <t-T> 


< 1 . 2 > 
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The pilot neuromotor dynamic* can be represented approxi- 
mately by an adjustable first order lag < " " >, generated 

n 

naturally in the modeling by including the control rate as 
the "input" to be chosen to minimize the cost function <see 
equ, <l.l>>.The augmented state equations then become 

x<t> « *.<t> + B Q pi<t> ■* u> 0 <t> , <1 , 3> 

with 


% * I 


I A 


B I 
0 I 


B 



I ! 
I 1 I 


w 0 <t> » 


I w<t> I 
I 0 I 


<1 .4> 


It can be shown in this case that the optimal control p, <=*u> 
is the linear feedback law 


p<<-> = G *.<t> <1 .5> 

where G » C G_ , G 1, or equivalently 

X u 

T u<t> + u<t> * L *<t> = u <t> <1.6> 

n c 

where ; L = , and ft<t> is the estimate of 

state x(t). The feedback gain G is 

C = -R" 1 B l K 0 

where K 0 satisfies the Riccati equation 


< 1 . 7 > 
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K o + K o *o + Q o " K o »0 R_1 B 0 K 0 - 0 < 18 > 

in which 



! c t qc 0 


I 

I 

I 


* 

The state estimator consists of a Kalman filter and 
predictor/ and the effect of the motor noise v u <t> is 
included. We define as before the augmented state %<t) ®col 
Cx(t>,u(t)3 which satisfies 


*x-<t> » A, x<t) + B, u <t> + ui 1 <t> 
l 1 c l 


<1 .?> 


v u Ct> 1 

in which w«<t> - col C w<t> / “““ 3> B 4 « col C 0 / r~* 3 / 

I T 1 T 


n 


Tl 


and 


I 

A „ 1 A B 

I O -1 


I 


n 


The Kalman filter 
^,<t~T> from: 


W. 


W 

0 


I 


u 

r 2 

n 


generates the delayed 


< 1 . 10> 


estimated state 


where 


&<t~r> = A, *<t-T> + S.cTv 1 C y <t > 
i 1 l y p 

- C^t-r) 3 + B 1 u c <t-r> < 1 . 1 1 > 

CC.,03, and the error covariance matrix satis- 


fies 
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The predictor then generates tft<t> according to; 


a,12> 


V 


*<t> * £<t> +• e * C*t<t“T> - £<t-T>3 
><t> « A t ^<t> -f B t u c <t> 




Thus/ the human operator model remains linear, Based on 

* 

man-machine experimental results# each white observation 

noise v and white motor noise v were found to have a 
y i u 

o 

covariance proportional to the mean squared values of EC y^ 
2 

3 and EC u c 3 respect ively# i , e , # 


V, = p, it EC y? 3 i~l # 2# . . .m 

y 4 *■ * 


V S p IT EC u ] 
u c 


Therefore# with V is normalized with respect to EC y 


< 1 « 14> 
2 , 


with p^ **,01# it has a positive frequency power density 

2 

level of -20 db. When is normalized with respect to EC u c 
3 with P 0 -.003# it is approximately -25 db. Both of these 
values have been found to model a variety of simple tracking 
tasks for a sveral different plants (controlled element 
dynamics >. 
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Past approaches 
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In the optimal control model of the pilot, therefore* 

there are some pilot-related model parameters that need to 

be selected; i. e. * weighting of the cost function Q and R, 

time delay r * measurement noise v and motor noise v etc.. 

y 


There 
imental ly 
for a comp 
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weighting 
other mod 
motor nois 


have been some frequency domain 
identifying these parameters d 
ensatory task. In this task only 
yed to the human controller and 
of error measurement is require 
el parameters* time delay • mea 
e* all have been identified (Ref 


methods' of exper- 
eveloped primarily 
the error signal 
therefore only the 
d. Combined with 
surement noise and 
. 1 , 6 * 0 * 19 * 17 ). 
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Now we want to identify the optimal control model 
parameters for pilot performing more complex task.** and 
select the weightings of error# error rate# and any other 
displayed variables appropriate in the task. We may identify 
the model parameters using bath frequency domain and time 
domain techniques. The former method is a classical 
approach# and the later is a new one developed by D. K. 
Schmidt# with emphasis on determining the cost function 
weightings (Ref. 14). We introduce these two methods as fol- 
lows. 

4 

3- A.* Time d oma in identification 
From section 1# we assume 

u » - G uVo <3.1> 

where G , G are the optimal control gain related to the 
X u 

weighting of cost function J, is the estimate of the state 
x, and v is the motor noise. Let € - x- ft , then 

U X 

° = V " G X 6 X + G u u " G u v u < 3a > 

It is assumed that the whole system is stationary and satis- 
fies the ergodic hypothesis. Taking data on u# u# and x from 
simulation over an appropriate time period# one can obtain 

1 G T 1 

! * ! = C M‘ - M 0 ]“* N. <3.3> 

G ' 1 * u 


where 



r 


M 4 * 


- 9 


[ E<xx T > E<xu T > j 
I E<ux T > E<uu T > I 


original PAge is 

F P0 °f> QUALITY 


M ss 

2 


E<xe3> ECXV, 1 ") ! 

A v I 

E<ug^> E<uv^> J 


N. « 

u 


I ' T ' 

I E<xu 1 > I 

I • t * 

I E<uu’> i 


<3 . 4> 


Here we rote that the estimation error e and motor noise v 
are unmeasurable. However we can still proceed with 
knowledge of the pilot model structure combined with a rea- 
sonable assumption on pilot time delay and the covariance of 
measurement noise and motor noises. 

In fact* we can obtain the following relations: 


E<xv y > ~> O 

E<U<£> x o 


and approximately u.e have 


<3 . 5> 


_ . 1/2 .. . 005ir 2. 

E< °V S T = E<U C > 

n n 


<3. 6> 


Her© we assume V has a normalized value of ~20db. E<x e > 

U X 

can be determined using the properties of the augmented 
state Kalman filter and a least mean square predictor. That 
is, since E<&e >=0, 


E<xeJ> = E< 6)< J> 


<3 . 7> 

found. We can then find that <Ref. 7> , recalling that C v6 


3=col Cx,u3 
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_ A,r aTt t Ajt Alt 
E<e e'> - e 1 E.e 1 + fe 1 W,e 1 dt <3.S> 

Tv A* i q I 

00 At ^i T T —1 A-t rTi 

+ X e e 1 SjCTv l C 1 S 1 e 1 e H X dt 
O i l y i i - 

Her©/ the A^/Ej/Wj have been defined in ©qu. < 1 . 10>-<1 . 12> , 

and are also functions of r <=“•> , and 

n b 

u 


I 

A » * A 

I zL 

I T 

. n 


B 

zi 

T 


<3 . 9> 


W© can now us© an iterative method in the identifica- 

j 

tion procedure. Try an approximates to calculate E< e e > 

n X- 

used to fit the gain vector C G ,6 3, then iterate the 

xu 

procedure until the approximated approaches its fitted 
va 1 ue . 


The final step is to find the meaningful values of Q 
and R which corresponds to the experimental ly determined 
optimal gains CG , G, 3 (see equ. <1.7> and <1.8>). In gen- 
eral the inverse solution is not unique. So in our identif- 
ication procedure me can arbitrarily choose the weighting on 
control rate <R> equal to 1, and then use a Quasi -Newt on 
method Csee appendix* to identify the other weightings in Q 
to minimize the modeling error function: 


N G X i % G -8 2 

j s < — i — h)* + <-y — 

1 i-1 *G °G 

x . u 


<3. 10> 


where G and G are estimated mean values/ from several of 

x 


i- 
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simulation runs; , or c are the experimental standard 


deviations corresponding to each gain; 6^ , C y are the gains 
from exercising the model corresponding to some value of Q. 


Natural ly, this method emphasizes the weightings in the 
cost function, because we select the time delay, measurement 
noise, motor noise empirically. We can also select these 
parameters from the results of the next method, and then 
compare the results in both approaches. 

3 . 2 . Frequence! domain identification 

' First we derive in the frequency domain some charac- 
teristic equations from the theoretical pilot optimal con- 
trol model. From equ. < 1 . 1 1 >-< 1 . 13> , we can obtain 

u <s> = L ■ &<s> = H<s> C Y < s > +V <s> 3 <3.11> 

c e y 

where 

L = CL, 03 
e 

H<s> 

with 

A = Aj - S 1 CjVy 1 C 1 <3. 13> 

and A as in equ. <3.9>. Therefore we can consider the pilot 
model block as Fig. 3-1. 

Considering now the pursuit task with single input, in 


T <sI-A >t 

= -e t L C <sI-A> /e dt <sI-A> 

0 

+ si-A + b 4 l i" 1 s.cTy -1 

.1 e 1 1 y 


<3. 12> 




! 

L 

h 

; j! 

I i 

I \ 


| 

l 



i : 


ij 
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which the subject observes the input command as well as the 
system output* < and also their rates implicitly). For exam- 
ple» we take the error e(t)# error rate eit) # input command 
i(t> and command rate i(t) as our observations (measure- 
ments). Then one can consider a closed-loop system # as 
shown in Fig. 3-2. Referring to the block diagram of Fig. 
3-2# we have 

E(s) - I(o) - M(s> 

. M( 5 ) a U(s>- F(s) ..(3, 14V- 

and 

U a <E+N 1 > + <sE+N a > H 2 + 

n 

<I+N 3 > H 3 + <sI+N 4 > H 4 + <3 . 15> 

where N« , N., and N are the measurement noise and motor 

1 t£ o u 

noises. Here obmiting the symbol <s>, we can derive 


r 





Figure 3-1 Equal pilot OCM model block diagram 
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E = 


F<H ~+sH.> 

I r I - - 

u ’ t s+l J 

n 


7^1 CH lVWWWV 

15 


1 . + 


F<H 1 +sH a > 

T S+l 
Tl 


<3. 16) 


H ^ +3H P +H«+sH 4 1 

I C- ^ a 3 + T '- ll ' T CH.N.+HAm+H.N.+N] 


U « 


T S+l 
Tl 


r s+l 11 2 2 ' 4 4 "W 

T> 


1 . + 


FCH^+sHg) 
r s+l 

Tl 


(3. 17) 


Now assume that the system input i(t)» measurement 
noise n« , n~.» n 0 ) , n . t and motor noise n. are uncorrelated with 
one another. Then we have the following power spectrum rela- 
tions: 




H j +sH a +Hg+sH^ 

T^S+I 

I-/U , r\ r <1* . . (o) 

F<H« +sl-U> ix 


1 + 


r s+l 
n 


<3. 18) 


$ < Q> = 

uu 


<t>. <o> = 

le 


1 - 


1 + 


F<H 3 +sH 4 > 

T n s+ 1 


F<H 1 +sH p > ii 

T S+l 

n 


<l> . . <o> 


Hj+sHg+Hg+sH^^ 


t s+l 
*n 


i + 


F<H 1 +sH 2 > 


T S+l 

n 


♦ + 


1 . 


T S+l 
5 


<3. 19) 


1 + 


F<H l +sH 2 > 


T S+l 

n 


! h i !%«/“> + ! h 2 ! 2 v 2 <^ + 


I 


!kJ 2 * „ <o) + !h.! 2 ^ ^ <o>+4> <o>l <3. 20) 

I 31 n 3 n 3 I 41 n 4 n 4 n^ , 

Now we define the equivalent describing function Y <j©> as 

D 

*eq 
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„ *lu<»> 


r P < jo > - 

'eq , ie 


<3.2i> 


Combining equ. (3.10) and <3. 19)» uie hav# 


. € jo> * 


*j U «» 

1® 


Hi+sH a +Hg+sH 4 

r s+1. 
n 

F<H 3 +sH 4 > 

1 T S+1 

n 


< 3 . 22 ) 


which is the same result as Hess / (Ref. B>* with observation 
of error only in compensatory task <i.e., Hg^H^O >. 

We may also define the controller remnant-correlated 

power spectrum d> uo as the part of the pilot input power 

r 

spectrum induced by the remnant < measurement noise and motor 
noise) , i , e . , 


♦uu <0> ’ *uu <8> - 
r 


H ^ +sHg+Hg+sH 4 

t s+1 

33 

F<H 1 +sH_> 

1 + 

T S+1 
n 


^ <o> 


T s+1 
n 

F<H 1 +sH 2 > 

n 


I I is 

| , ♦ <o> 

I II n^n^ 


+! h p!\ „ <®>+[h ' a <* „ <©> 

I 2 1 n^n^ ® I ^3^3 

+! h J 2 * ' <©>+♦_ „ <»>! 

1 41 1*4*4 Vu I 


<3 . 23) 


Equ. <3 ..22) and <3.23) are the character istic equations for 
a single input pursuit task. Using the same procedure, it 


can be extended to other complex tasks. 


In our certain selection of input command (see section 


i 


- 16 - 


ORIGINAL PAGE IS 
°P POOR QUALITY 


4>, the input power spectrum ♦ * s z**' 0 at other than 
input frequencies. So analyzing the simulation data* uj* can 
calculate the equivalent describing function only at the 
input frequencies! and obtain the controller remnant- 
correlated power spectrum at other than input 
frequenc iesi i. e. i 

« i nput frequenc i es 


v ✓ . s 

Y p c _ < ' i °k > “ 


eq 


ie k 


^uu < " < * > i > ** ^uu^l^ input frequencies <3.24> 

r 

For computing the power spectrum, we calculate the sig- 

i 

nals correlation first, and then fast fourier transform it 
to obtain the power spectrum. 


We can identify these model paramet *rs using theoreti- 
cal modeling to match with the simulation result. For fit- 
ting these parameters , we can also include some other 

important effects which can be obtained during the simula- 

2 

tion, such as the variance of pursuit error E< e > and the 

2 

variance of controller E< u >, Finally we also use Quasi- 
Newton identification procedure to minimize the modeling 
error defined below: 


N -*■ N ^ 

i 1 G r G i 2 i 2 p r p i 2 

J« - rp 2 <~z — L > + rr- E — L > 2 

N l i=i *g. N 2 i=l a P. 


N ^ N a 

1 3 R i~ R i 2 1 4 S i” 3 i 2 

+ rr 2 C- 1 + rp E C— -> d 

N 3 i«l a R. N 4 i=l a S. 


<3 . 25> 


where 


t 
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N. «no. of valid measurement it) the i group 
C; “magnitude of the i* n describing function point 
to be matched/ db. .. 

P. “phase shift of the * describing function point 
to be matched/ deg. 

R. “controller remnant-correlated power spectrum, 
of the.i vh frequency point to be matched/ db. 

S. » i tr> variance score to be matched 
“standard deviation of experimental data 

" A " : indicates model prediction 
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4. 1,. Introduction 

The experiment with the pilot in the loop is done in 
the Flight Simulation Lab.! with the use of the minicomputer 
i CRT 1 and control stick etc.. In our casei we display the- 
input command i(t) and the system output m<t) on the CRT 
screen! the human pilot observes the display and determines 
his command into the control stick! in his attempt to null 
the error between input command i(t) and system output m(t>. 
Finally the analog signal from the control stick is con- 
verted to digital and is input to the dynamic system being 
simulated numerically by the minicomputer. This closed-loop 


system 

with human in th 

e loo 

p is 

shown in Fig. 

4 

-li which 

may 

b 

e compared to the 

OCM 

mode 

1 already shown 

in 

Fig. 

3-2. 

The 

ta 

sk to be treated 

here 

is 

the 

pursuit tr 

ac 

k ing 

task 

with 


a single input. 

The 

CRT 

and 

control stic 

k 

in the 

lab 

are 

sh 

own in Fig. 4-2. 

Fig. 

4-3 

shows the CRT d 

is 

play 

f or- 

mat# 


in which the distance 

between 

line A and 1 

in 

e C is 

the 

input 

command i(t>* the 

line 

s B 

and 

C represent 

sy 

stem 

out- 

put 

m 

(t)# and line C i 

s the 

zer 

0 re 

iference line 


During 

our 


experiment# we simulate two simple systems and 








Figure 4-2 Simulation equipment; CRT and control stick 
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respectively. 




4. g,* Selection of the sinusoidal input 


In the simulation# we use a sinusoidal random-appearing 
input* mechanized as the sum of sine waves, i. e. » 

N 

i<t> « E A.. sin<o L ,t -t* > <4,0 

k=i K K K 

Here the are chosen to be nbn-commensurable <no frequency 
is an integral multiple of another), and roughly evenly 
spaced on a logrithmic scale, In addition the are 

selected so that in a finite experimental run length, all 
the constituent sine waves in i(t) will have completed an 
integral number of cycles. Finally the are chosen to lie 

within the range of interest of human response workj i. e. # 

0.1 < o. < 20. <4 , 2> 

~ k ~ sec 

The autocorrelation of the sinusoidal input is 


2 

♦ii <T > = S~ COS< °k T > 
and its covariance is 


<4. 3) 


a?. * S ~ <4 . 4> 

1 k = l 

Sc when we use a sinusoidal input as a random-appear ing 
input, we should appropriately select these amplitudes A^ to 
match the frequency distribution of power in the real random 
input power spectrum being approximated. 
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For example* we want to use a sinusoidal input to 
replace- a real random input command 8 c , which satisfies: 


1 ; 1 * 0 I 

I c II 0 1, I I B c I I 0 I 

II l I -2.23 *-3 111 1111 

I 0 c I I 0 c I 


<4 . 5 > 


Where w is a white noise A 'N<0* 13, 3), Taking the Laplace 
transform: 


%* 8 <s> + 3s 0 <s> + 2.25 8 <s> » w<s> <4.6> 

c c c 

or 


then 


$<JZl 

w<s> 


la. 


<S+1 .5>‘ 


<4 . 7 > 


l l 

! i. ! . 13.5 

<C> « I 2 1 ♦ U uu < ® > 3 2 2 

c 0 c I < Jo+l , 5> 2 1 UU> +2.25> a 


<4,8> 


To select the sinusoidal input to match the frequency 
distribution of power with the real random signal* first we 
define the fraction of power of real random signal 8 <t> 
<Ref, 13> between 0 ( « as 


o 




p* <«>do h* <o>d» 

F e • <“> - ;« 

C C if ^ ^ . & ■ 


<4 . 9 > 


”1 *0 c B c (a>da 


8 8 
c c 


And we also define the fraction of power of the sinusoidal 
input i(t> as 




n<o <o> 

1 2 

1 5 ? 6 “ 

2 * , H k 

k”l 

1 N g 

i E , A k 

k“l 


n<o <©> 

1 71 2 

2 E A k 
a k=i k 


<4. io> 


or 


ii 


I 
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Certainly we select a . Final ly, we select the 

ii 0 c e c 

appropriate A^ to match the fraction of power spectrum of 
sinusoidal input with that of the real random input. It is 

shown in Fig. 4-4. Table 4-1 shows the parameters selected. 

* 

During our system dynamic computation, we choose an 
integration time interval dt equal to .05 sec and. the number 
of data samples equal to 1024< 2 10 >, for the sake of con- 
venience for the fast Fourier transform. The period time T 
during which data is taken is equal to 51.2 sec, which 
satisfies the following relation 

T “ N At <4 . 1 1 > 

and define 

= = .1227185 3C “ < 4 . 12 > 

' 0 T sec 

We then choose some integers n^ which are non-commensurable 
to obtain the input frequencies, or 

w k * "k <4.13> 

and also are to be equally spaced over the logarithmic 

scale and satisfies the requirement in equ. <4.2>. 

4.3. Simulation result 

In our experiment, the sinusoidal input shown in Table 

4-1 is used to represent a true random input. We simulate 

k k 

two simple plants — and ~ in both a pursuit and a compensa- 

s 

tory task. The results for the pursuit task are shown in 
Fig. 4-S, 4-6 and Table 4-2 respectively. Fig. 4-7, 4-8 show 


000*1 
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Igure 4-4 Comparison of frequency distribution of power In sum of sine waves with random signal 
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Table 4-1 Sum of sinusoidal Input command 


lterm 

k/s dynamics 

k/s* dynamics 

mean value 

> 

deviation 

mean value 


squere of 
error 

.26 

.06 

.61 

.13 

square of 
error rate 

4.06 

.74 • 

5.08 

1.03 

I 

EHa 

3.05i 

.72 

26.78 . 

8.65 


Table 4-2 Measured human performance' 
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Figure 4-5 Measured human pilot equivalent describing function 
and controller remnant-correlated power spectrum 
for k/s dynamics 
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<jJ Y&d/sec 


{ H l M U 


Figure 4-6 Measured human pilot equivalent describing function 
and controller remnant-correlated power spectrum 
for k/s* dynamics 
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Figure 4-7 Comparison of compensatory task with other results 
for k/s dynamics 
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Figure 4-8 Comparison of compensatory task with other results 
for k/s 2 dynamics 
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the result* tor the compensatory task* compared with other 
result* in the literature (Ref. 1<21»22). 

From these later result* we conclude that our qxperi- 

mental technique and software are correct and that we may 

* 

proceed to more complex task*. 

4. Discussion of Experimental Da ta 

A fixed time interval (.05 sec) and a second order 

Runge-Kutta integration are used in the real time Simula- 

tion. The whole discrete closed-loop system is sampled at 20 

&W£ , -lfi , ^=125 > which is much larger than our maximum 

sec sec * 

» p 

input frequency <=17 >. So no aliasing problem is 

expected for this sampling rate, but at slower rates (say 

less than 10 it must be considered, 

sec 

Our simulation results shows that there i* a little 
difference between pursuit and compensatory tasks. The mag- 
nitude of the pilot describing function is smoother over the 
whole frequency rang# in the pursuit task than in the com- 
pensatory task. We also have phase lead in the range of low 
frequencies in the pursuit task, but not in the compensatory 
task. Also, a lower error variance occur* in the pursuit 
task. since the subject has more observations, bath input 
and system feedback, and therefore the task is "easier" than 
the compensatory task. 
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We want to identify the pilot dynamic* and task objec- 
tives (weightings) from simulation data similar to that 
described* using both time domain and frequency domain 
methods# and compare the results, 

We propose using the Quasi-Newton identification pro- 
cedure (see appendix) to identify these model parameters* 
such as weighting of cost function* time delay# measurement 
noise and motor noise etc. . 

Our proposed technique can be extended to other more 
complex tasks* for example# landing approach. During that 
task* we have more observed variables such as altitude# 
attitude* angle attack and velocity etc. > and more controll- 
ers such as thrust and elevator etc. # than those we have in 
our pursuit task here. Our proposed technique is still 
available# but more weightibgs of the cost function must be 
selected# due to the increasing observed variables and con- 


trol iers. 
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APPENDIX 

A. Quasi-Newton identification procedure (Ref, 15# 16» 17) 

A. 1 Minimization scheme 

The QN scheme is generally implemented to minimize the 
following scalar modeling error: 

N 2 

J « E M.ef <A~1> 

i=l x x 

where is the difference between the i*^ 1 measured data 
point and the corresponding model prediction, w^ is a weght- 
ing coefficient. Or in matrix form: 

J - e T We <A-2> 

with e^col C e^ e g W=diag [ ], 

For a trial set of model parameters p^ , we have its 
corresponding modeling error 

™ e i^ e i <A~3> 

For a new set of parameters p =p,+Ap , we obtain a new 

d X 

modeling error 

$ 

J 0 - <©j+Ae> T W Ce^+Ae) <A-*4> 

— e^We^+2e^NAe+Ae^"h4Ae 

Using perturbation theory ( we can get approximately linear 

perturbed equations in the model parameters. Thus 


Ae ~ Q Ap 


<A-5> 


where q<i,j>-g — can be obtained by the method of perturba- 

P j 
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tlon. Now th» modeling error can be expressed as 

J 2 * J 1 +2eJwQAp-fAp T Q T WQAp <A-6> 
Next me want to find the parameter vector change required Ap 
for minimizing the Jg , There 

t 

S'j 2 ->minimum " 0 = M^WhoAp <**> 
Thusi the following change in the parameter vector yields 
minimum modeling error/ given the initial vector e^ and the 
assumption of linearity* 

Ap » - CQ T WQ3" 1 Q T Ue 1 <A-8> 
A. 2 Sensitivity analysis 

In addition to obtaining the best match to a given set 
of data/ we may also wish to determine some measure of the 
reliability of the identified parameter values. A qualita- 
tive indication of parameter estimation reliability can 
often be obtained through sensitivity analysis relating 
changes in the scalar matching error to perturbations in the 
model parameters. In general* estimates of parameters that 
have a high impact on the modeling error can be considered 
more reliable than estimates of parameters having a smaller 
impact. 

If model predictions are linear in the parameters* as 
assumed in the foregoing treatment* we may analytically 
derive the sensitivity of the scalar modeling error to per- 


turbations in model parameters about the optimal (best 


I 
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matching) set. One may compute the sensitivity to a given 
parameter with the remaining model parameters held fixed* or 
u/ith remaining parameters reoptimized. The latter measure 
provides a more accurate reliability measure because it 
accounts for the potential tradeoffs that may exist among 
parameters in terms of matching the data. 

Let e^ be the modeling error when parameter set p is 
optimized at p^ , then: 

Ap = -rQ T WQ3~ 1 Q T We <> = O <A-9> 
Next let us assume that the incremental error arises from a 
non-optima l choice of one single parameter p^ . With the 
remaining parameters fixed at their optimal values* the 
resulting incremental error is 

Ae = Ap^ <A-10> 
where q^i col. of Q. We define the subscript "r" to 
indicate vectors and matrices that remain when rows and 
columns corresponding to the i*^ model parameter are 
removed. The expressions of re-optimizing the remaining 
model parameters can be obtained: 

Ap r = -CQ^WQ r 3“ 1 Q^Wq.Ap. <A~11> 
Comparison of the elements of the vector Ap^ with p A reveals 
the joint tradeoff between p i and the remaining model param- 
eters. 

To compute the effect on the modeling error J of a 
change in p 4 , with the remaining parameters re-optimized. 
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we construct a new vector Ap which is the composite of 

p. and p . This vector is defined as 
l r 

Ap “ V Ap^ .<A-12> 

where V is a column vector that has a value of unity for ths 1 
i element and values for remaining elements as determined 
from Ap . Then we can obtain the corresponding modeling 

r » 

error for re-optimizing p^ with the change of one single 
parameter Ap^ 

AJ - J-J Q - 2eJwQUAp.+V T Q T WQ0<Ap i > 2 

- 0 T Q T WQ0<Ap i > 2 <A-13> 

T 

The term 2egWQVAp i is zero, because e Q is corresponding to 
the optimized modeling error. Therefore the change in model- 
ing error varies as the square of the change in the parame- 
ter value. Hence we can obtain the sensitivity for each 


parameter. 
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